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o 

The electronic structure of the ferroelectric crystal, NaN02, is studied by means of first- 
principles, local density calculations. Our ab-initio, non-relativistic calculations employed a 
local density functional approximation (LDA) potential and the linear combination of atomic 
orbitals (LCAO). Following the Bagayoko, Zhao, Williams, method, as enhanced by Ekuma and 
Franklin (BZW-EF), we solved self-consistently both the Kohn-Sham equation and the equation 
" q giving the ground state charge density in terms of the wave functions of the occupied states. We 

found an indirect band gap of 2.83 eV, from W to R. Our calculated direct gaps are 2.90, 2.98, 
i— h 3.02, 3.22, and 3.51 eV at R, W, X, T, and T, respectively. The band structure and density of 

+3 states show high localization, typical of a molecular solid. The partial density of states shows 

that the valence bands are formed only by complex anionic states. These results are in excellent 
agreement with experiment. So are the calculated densities of states. Our calculated electron 
effective masses of 1.18, 0.63, and 0.73 m Q in the T-X, T-R, and T-W directions, respectively, 
show the highly anisotropic nature of this material. 

Pacs Numbers: 77.84.-s, 71.20.-b, 71.20.Mq, 71.20.Nr 

O 
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I. INTRODUCTION AND MOTIVATION 

(N 
> 

Sodium nitrite is one of the ABO2 group with several interesting properties such as ferroelec- 
y— \ tricity, paraelectricity and piezoelectricity. The ground state of NaN02 has one of the simplest 

f — structures in the ABO2 group. Since the discovery of its ferroelectric property in 1958, by 

Sawada et al. [1], there has been increased interest in diverse properties of NaN02, ranging from 
structural changes to dielectric, electronic, thermal, electrical and elastic properties, as well as 
its non- linear optical spectra [2J. Despite this noted interest in ferroelectric NaN02, there has 
been very little study of the electronic band structure of the crystal [3] . 

One of the earliest electronic band gap measurements of ferroelectric NaN02 was by Asao 
J> et al. [4] who found band gaps of 2.3, 2.6, 3.1, and 3.7 eV for different current directions in 

their resistivity work. Also, using the same method, Takagi and Gesi [5] found a band gap of 
r> 2.4 eV. These authors did not specify whether the gaps in question are direct or indirect. The 

infrared absorption spectra study of Sidman [6] at a temperature of 77 K gave a band gap of 
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3.22 eV while that of Verkhovskaya and Sonin [7], at 293.15 K, led to a value of 3.14 eV. The 
semi-empirical LCAO work of El-Dib and Hassan [8], for room temperature (293.15 K), found 
a band gap of 3.24 eV. The scanning electron microscopy, diffuse reflectance spectroscopy, and 
electrical measurements study of Balabinskaya et al. [9] , at a temperature of 300 K , reported 
an absorption edge of 3.0 eV. 

One of the earliest theoretical studies of the band structure of ferroelectric NaN0 2 was the 
Xa (a = 0.75) exchange calculation of Kam et al. [10] who found a direct band gap of 3.45 
eV at r and an indirect band gap of 2.0 eV. The Full-potential linear muffin-tin orbital (FP- 
LMTO) method of Ravindran et al. [2] found an indirect band gap of 2.2 eV. The generalized 
gradient approximation (GGA) work of Wang et al. [11] utilized a full-potential linearized 
augmented plane wave (FL-APW) method to calculate a band gap of less than 2.1 eV. The 
LDA approach of Zhuravlev and Korabel'nikov [12], using the nonlocal Troullier-Martins (TM) 
pseudopotentials and the Slater exchange potential with a correlation correction, led to a direct 
band gap of 2.50 eV at W. Zhuravlev and Poplavnoi [13], also using nonlocal Troullier-Martins 
(TM) pseudopotentials and the Slater exchange (a = 1) potential with a correlation correction, 
found a minimum forbidden band gap of 3.07 eV. The LDA result of Jiang et al. [14], using 
an orthogonalized linear combination of atomic orbitals (OLCAO) method, found an indirect 
band gap of 2.95 eV from S to R symmetry points. The periodic Hatree-Fock (PHF) method 
of McCarthy [15] overestimated the gap by a factor of 2-3. Henkel et al. [16] employed their 
plane-wave and Gaussian basis set and an Xa exchange and reported a direct band gap of 2.70 
eV. These results suffice to see that theoretical calculations have not resolved the under- or 
overestimation of the band gap of ferroelectric NaN0 2 - 

In light of the preceding overview, the aim of this study is to attempt to obtain the mea- 
sured, fundamental band gap and other electronic properties of ferro-NaN02 with ab-initio, 
self-consistent LDA calculations. The confirmation of our predictions of the band gaps and 
other properties for BaTi0 3 [17], Ti0 2 [18], CdS [19], c-Si 3 N 4 [20], c-InN [21], and SrTi0 3 [22] 
is a basis for the above presumption. Further, the mathematical rigor of the method [18-24] 
suffices to expect it to lead to much better results. 



II. METHOD OF CALCULATION 

Our calculations employed the Ceperley and Alder [25] local density functional potential as 
parameterized by Vosko, Wilk, and Nusair [26]. We implemented the linear combination of 
atomic orbitals (LCAO), using Gaussian functions for the radial parts. We utilized a program 
package developed at the Ames Laboratory of the US Department of Energy (DOE), in Ames, 
Iowa [27]. Our calculations are non-relativistic and were performed using low temperature (120 
K) lattice constants. The distinctive feature of our approach resides in the implementation of the 
Bagayoko, Zhao, and Williams (BZW) method, as enhanced by Ekuma and Franklin (BZW-EF), 
consisting of concomitantly solving self-consistently two coupled equations as explained below. 
One of these equations is the Schrodinger type equation of Kohn and Sham [28] , referred to as 
the Kohn-Sham (KS) equation. The second equation, which can be thought of as a constraint on 
the KS equation, is the one giving the ground state charge density in terms of the wave functions 
of the occupied states. 
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The essentials of the method follow. Beginning with a minimum or small basis set capable of 
accounting for all the electrons in the system under study, one performs ab-initio, self-consistent 
calculations. This minimum or small basis set is constructed using orbitals resulting from self- 
consistent calculations for the atomic or ionic species in the solid. As noted below, these species 
are Na 1+ , N 1+ , and O 2 " for NaNC>2. Subsequently, this basis set is augmented with one orbital 
and self-consistent calculations are done. This new orbital and others that may be added later 
are outputs of the self-consistent studies of the ionic species. The occupied energies of the two 
calculations for the solid are compared numerically and graphically. These occupied energies 
from Calculations I and II are generally different. A third calculation is carried out after adding 
another "ionic" orbital to the basis set. We note that adding an orbital means augmenting the 
size of the basis set by 2, 6, 10, or 14, depending on the s, p, d, or / character of it, respectively. 
Again, the occupied energies from Calculations II and III are compared. This process continues 
until the occupied energies from a calculation, i.e., N, are found to be identical to those from 
Calculation (N+l) immediately following it, within our computational uncertainties of 50 meV or 
less. At that point, the calculations are completed and the results from Calculation N represent 
the physical description of the system under study. The basis set for Calculation N is referred 
to as the optimal basis set. Calculation (N+l) and others, with larger basis sets that contain 
the optimal one, provided linear dependency is avoided, reproduce the occupied energies from 
Calculation N. However, by virtue of the Rayleigh theorem [19,23], these calculations produce 
some unoccupied eigenvalues that are lower than those from Calculation N, due to a non-trivial 
basis set and variational effect [23]. 

All the self-consistent calculations known to us carry out the iterative procedure that involves 
both the Kohn-Sham equation and the equation giving the ground state charge density. However, 
the methodical increase of the size of the basis set, as done in the BZW-EF method, entails 
changes in the radial and angular characteristics of the basis set as well as an increase in its size. 
It is in this sense that the BZW-EF method solves self-consistently the two equations in question 
to obtain, in a verified fashion, the minima of the occupied energies. The difference between 
our method and single trial basis set calculations may be best understood by noting that any 
serious radial, angular, or basis size deficiencies that may exist in a single trial basis set cannot 
be remedied by the iterative process noted above. In the case of the BZW-EF method, such 
deficiencies are corrected as the methodical augmentation of the basis set is carried out. The 
attainment of the minima of the occupied states signifies that there is no deficiency left to correct. 
When a larger basis set that contains the optimal one is utilized to carry out a calculation, the 
charge density, the potential, and the Hamiltonian resulting from self-consistency, as well as the 
occupied eigenvalues, are the same as the corresponding ones obtained with the optimal basis set. 
Hence, the Rayleigh theorem applies to the outputs of this calculations whose occupied energies 
are the same as those obtained with the optimal basis set. However, some unoccupied energies 
from the larger basis set are generally lower than (or equal) to their corresponding ones obtained 
with the optimal basis set. The lowering of these unoccupied energies cannot be ascribed to 
physical interactions, as the Hamiltonian did not change. It results from the non-trivial and 
well-defined basis set and variational effect. 

The enhancement of the original BZW method is in the methodical increment of the basis set. 
This enhancement leads to adding the polarization {p, d, or f) orbitals, for a given principal 
quantum number, before adding the spherically symmetric s orbital (see Table 1). These addi- 
tional unoccupied orbitals are needed to accommodate the reorganization of the electron cloud, 
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including possible polarization, in the crystal environment. For valence electrons in molecules 
to solids, polarization has primacy over spherical symmetry. 

While there are no specific rules governing the order of adding orbitals, our experience has 
been that orbitals often need to be added to the heaviest element before the others. Among 
these others, preliminary studies of charge transfer serve as a significant guide. In the case of 
NaNC>2, even though oxygen is a little heavier than nitrogen, the neon configuration oxygen 
approaches after gaining 2 electrons intimate more stability than that of N 1+ . The ultimate 
aim of the addition of orbitals (i.e., basis functions) is to allow an optimal reorganization of the 
electronic cloud in the solid environment as compared to atomic or ionic ones. 

Computational details germane to a replication of our work follow. Ferroelectric sodium 
nitrite (ferro-NaNC^) has an orthorhombic structure [29]. It is in the space group C|v _ IMM2, 
with space group number 44 and Patterson symmetry lmmm [30]. The ferro-NaNC>2 unit cell 
contains eight atoms: four (4) cations and four (4) anions whose positions are as indicated 
between parentheses: Na: (0,0.5881,0), (0.5,0.0881,0.5); N: (0,0.1224,0), (0.5,0.6224,0.5); and O: 
(0,0,0.1962), (0.5,0.5,0.6962), (0,0,0.8038), (0.5,0.5,0.3038). 

We carried out five (5) different calculations utilizing five (5) different basis sets in search of 
the optimal one. Table 2 contains the basis sets utilized in the five (5) self-consistent calculations. 
Methodical increases of the basis set led to calculation IV as the one with the optimal basis set; 
i.e., the calculation yielding the minima of all the occupied energies. Calculation V does not 
lower any occupied energies as compared to Calculation IV. Hence, the electronic structure and 
related properties presented here were obtained with basis set IV, the optimal basis set. 

Our self-consistent computations were performed utilizing experimental lattice constants of 
3.518 A, 5.535 A and 5.382 A for a, b, and c, respectively, as obtained at 120 K [31]. Neutral 
charge calculations (Na°, N° and O ) were carried out. Calculated charges were found to be +1 
for Na, +1 for N and -2 for O. Then, we performed self-consistent calculations for Na 1+ , N 1+ and 
O 2 " to get the input quantities for the calculations for NaN02- The radial parts of the atomic 
wave functions were expanded in terms of Gaussian functions, employing a set of even-tempered 
Gaussian exponents. 

The self-consistent atomic calculations provided trial atomic potentials for Na, N, and O, 
respectively. These potentials were used to construct the initial potential for ferro-NaN02- 
We used 16 Gaussian functions for the s and p states and 14 for the d states for Na and 
N, respectively. We utilized 17 Gaussian functions for the s and p states and 15 for the d 
states for O. A mesh of 48 k points, with proper weights in the irreducible Brillouin zone, 
was employed in the self-consistent calculations. A total of 121 weighted k-points were used in 
band structure calculations and a total of 147 weighted k-points were employed to generate the 
energy eigenvalues for the electronic density of state calculation. The computational error for 
the valence charge was about 0.073343 for 52 electrons, a little more than 10" 3 per electron. The 
self-consistent potentials converged to a difference around 10" 5 after about 60 iterations. 

III. CALCULATED ELECTRONIC STRUCTURE 

Our calculated, ab-initio, self-consistent bands for ferro-NaN0 2 (Calculation IV) are as shown 
in Fig. 1 and the comparison plot of the electron energy bands for basis set IV (solid lines) and 
basis set V (dashed lines) are as shown in Fig. 2. As is apparent in this graph, the eigenvalues 
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of the occupied states totally converged, within computational errors of about 0.050 eV, showing 
clearly that Calculation IV is the one with the optimal basis set. 
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Fig. 2. The calculated electronic energy bands of ferro-NaNC>2, from Calculations IV (full 
lines) and V (dashed lines). The occupied energies from Calculations IV and V are practically 
identical. 
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Fig. 3. The calculated density of states (DOS) of ferro-NaN02, obtained with the BZW optimal 
basis set (i.e. Calculation IV). 
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Fig. 4. The calculated partial density of states (pDOS) of ferro-NaNC-2 as obtained with the 
BZW optimal basis set (i.e. Calculation IV). 
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Table 1. Search for the optimal basis sets (Orbital added is in bold), as per the BZW method, 
for the description of the valence states of ferroelectric sodium nitrite (ferro-NaNC^)- The 
optimal basis set is that from Calculation IV. Na is Sodium, N is Nitrogen, and O is Oxygen. 
The maximum of the Valence band occurred at W, while the minimum of the conduction band 
occurred at R. NaN02 is an indirect band gap material. 



Basis 
Set 


Na, JN 
and O2 
Core 


Na Valence 


N 

Valence 


O2 Valenc 


Total No. 
:eof Valence 
orbitals 


Band Gap (eV) 














Direct at 

r 


Indirect 


I 


Is 


2s2p3s 


2s2p 


2s2p 


34 


2.40 


2.38 


II 


Is 


2s2p3s3p 


2s2p 


2s2p 


40 


2.70 


2.53 


III 


Is 


2s2p3s3p 


2s2p3s 


2s2p 


42 


3.18 


2.93 


IV 


Is 


2s2p3s3p3d 


2s2p3s 


2s2p 


52 


3.22 


2.83 


V 


Is 


2s2p3s3p3d 


2s2p3s 


2s2p3s 


56 


3.24 


2.83 



Table 2. The Calculated Eigenvalues (in eV) at the High Symmetry Points for ferro-NaN02- 
The eigenvalues are obtained by setting the energy at the top of the valence band, which occurred 
at W, equal to zero. 



X 


-25.11 


-25.11 


-21.03 


-21.03 


-21.01 


-21.01 


-20.97 


-20.97 


-20.05 


-20.05 


-9.24 


-9.24 


-6.73 


-6.73 


-6.21 


-6.21 


-6.08 


-6.08 


-2.51 


-2.51 


-1.80 


-1.80 


-0.01 


-0.01 


3.01 


3.01 


10.08 


10.08 


10.36 


10.36 


11.38 


11.38 


11.79 


11.79 


13.96 


13.96 


14.12 


14.12 


14.32 


14.32 


17.49 


17.49 


18.90 


18.90 


19.82 


19.82 


24.92 


24.92 


30.26 


30.26 














r 


-25.11 


-25.09 


-21.13 


-21.03 


-21.02 


-21.01 


-20.89 


-20.89 


-20.13 


-20.07 


-9.68 


-9.28 


-6.34 


-6.32 


-6.3 


-6.16 


-6.00 


-5.89 


-2.04 


-1.76 


-1.70 


-1.33 


-0.51 


-0.10 


3.12 


3.51 


6.66 


8.59 


9.31 


9.61 


10.31 


10.91 


12.53 


13.59 


13.86 


13.97 


13.97 


14.67 


14.78 


16.26 
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16.47 


18 


18.84 


20.17 


20.43 


20.86 


20.96 


21.97 


22.06 


31.43 














R 


-25.15 


-25.15 


-21.00 


-21.00 


-20.99 


-20.99 


-20.92 


-20.92 


-20.18 


-20.18 


-9.44 


-9.44 


-6.74 


-6.74 


-6.14 


-6.14 


-6.08 


-6.08 


-2.38 


-2.38 


-1.47 


-1.47 


-0.07 


-0.07 


2.83 


2.83 


9.18 


9.18 


10.12 


10.12 


11.02 


11.02 


12.02 


12.02 


13.4 


13.4 


14.46 


14.46 


15.62 


15.62 


15.73 


15.73 


17.96 


17.96 


22.35 


22.35 


24.19 


24.19 


24.99 


24.99 














W 


-25.15 


-25.15 


-21.02 


-21.02 


-21.00 


-21.00 


-20.98 


-20.98 


-20.06 


-20.06 


-9.28 


-9.28 


-6.84 


-6.84 


-6.13 


-6.13 


-6.08 


-6.08 


-2.41 


-2.41 


-1.79 


-1.79 








2.98 


2.98 


10.31 


10.31 


10.39 


10.39 


10.93 


10.93 


11.59 


11.59 


13.38 


13.38 


13.50 


13.50 


15.72 


15.72 


17.68 


17.68 


18.34 


18.34 


19.96 


19.96 


24.09 


24.09 


26.91 


26.91 














T 


-25.14 


-25.14 


-21.08 


-21.08 


-21.01 


-21.01 


-20.89 


-20.89 


-20.11 


-20.11 


-9.53 


-9.53 


-6.31 


-6.31 


-6.15 


-6.15 


-6.06 


-6.06 


-1.91 


-1.91 


-1.52 


-1.52 


-0.23 


-0.23 


3.28 


3.28 


7.33 


7.33 


9.74 


9.74 


10.95 


10.95 


12.71 


12.71 


13.11 


13.11 


14.43 


14.43 


16.59 


16.59 


17.99 


17.99 


18.98 


18.98 


19.60 


19.60 


20.71 


20.71 


25.00 


25.00 















Table 2 shows the energies at the high symmetry points. Figs. 3 and 4 show the calculated, 
total (DOS) and partial (pDOS) density of states derived from the bands in Fig. 1. The 
calculated peaks in the valence band DOS closest to the Fermi energy are at -0.25 ± 0.1 eV 
and -2.4 ± 0.2 eV. We found a third peak at -5.8 ± 0.1 eV which agrees rather well with the 
experimental value of -5.8 eV reported by Kamada et al. [32]. In the conduction band, we 
found relatively broad peaks whose centers are located at 8.25, 9.50, 10.50, 11.25, and 12.00 
eV, respectively. The widths of these broad peaks are about 0.5 eV. Caution is required in 
using the conduction band density of states at energies above 6 eV. As explained by Jin et al. 
[33] and Bagayoko et al. [24, 34], BZW-EF calculations are concerned with correctly solving 
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equations describing the ground states. Hence, while the low-laying excited states from these 
calculations have been found to agree with measurement, upper excited states are not expected 
to be correctly described. Specifically, the calculated optical properties of wurtzite InN [33] agree 
with experiment up to energies of 5.5 to 6.0 eV. The comparison of the bands from Calculations 
IV and V in Fig. 2 suggests that our calculated, excited states for NaNC>2 could be meaningful 
up to 20 eV. This situation could be a feature of molecular solids. 

Our calculated electron effective masses in the T-X, T-R, and T-W directions are 1.22, 0.64, and 
0.74 m Q , respectively. These values clearly corroborate the anisotropic nature of this material. 

Our calculations were performed at lattice constants for which experimental data are available 
for comparison purposes. However, the question arises as to what results could be obtained with 
a fully optimized crystal structure (i.e., equilibrium lattice structure). The optimization led 
to equilibrium lattice parameters of 3.529 A, 5.560 A and 5.391 A for a, b and c, respectively. 
The difference between our former (120 K experimental lattice parameter) and the latter (with 
optimized lattice parameters) eigenvalues is approximately 0.02 eV, with the energies obtained 
with the experimental lattice constants being lower. The two results are practically the same as 
their difference is smaller than our computational uncertainty of 50 meV (0.05 eV). 

IV. DISCUSSION 

An added dimension of the validity of the work reported here stems from the rigor of the BZW- 
EF method. Our calculations solved self-consistently the system of equations defining LDA. They 
are totally ab-initio and do not entail the use of two different theories; the sophisticated work of 
Zakharov et al. [35] utilized LDA wave functions as input in their quasi-particle calculations. 

Unlike single trial basis set (STBS) DFT and other calculations, LDA BZW-EF results agree 
with experiment and point to much more physical content of LDA eigenvalues than previously 
thought. As explained by Zhao et al. [23], STBS calculations do not avoid a well-defined 11 basis 
set and variational effect that has plagued electronic structure calculations since their inception, 
with emphasis on DFT calculations. 

Indeed, BZW-EF calculations begin with the minimum basis set that is just large enough to 
account for all the electrons in the system, i.e., NaN02 for this work. It then augments the basis 
set methodically [22-24] and carries out successive, self-consistent calculations with increasingly 
larger basis sets. Except for the first calculation, the occupied eigenvalues of a calculation 
are compared numerically and graphically to those of the calculation immediately preceding 
it. Ultimately, the occupied eigenvalues of two consecutive calculations, say N and (N+l), are 
found to be identical-within computational uncertainties of about 50 meV. The results from 
calculation N provide the physical description of the system and the corresponding basis set is 
the optimal basis set. While Calculation (N+l) yields the same occupied energies as Calculation 
N, some of the unoccupied energies from it are generally lower than corresponding ones from 
Calculation N. This extra-lowering of the unoccupied eigenvalues, for calculations with basis sets 
larger than the optimal one, is due to a basis set and variational effect. It is a mathematical 
artifact stemming from the Rayleigh theorem. For basis sets larger than the optimal one, the 
resulting charge density and potential are identical to those from Calculation N [23]. So, the 
above extra-lowering is not due to a physical interaction. 

From the density of states plots, our calculated peaks of -0.25 ± 0.1 eV and -2.4 ± 0.2 eV in 
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the valence band DOS agree with the x-ray photoelectron results of Calabrese and Heyes, [36] 
and the theoretical calculations of Ravindran et al. [3] and Zhuravlev and Korabel'nikov [12]. 
The observed sharp peak at -0.25 eV arises mainly from the hybridization between the O 2p , 
N 2p, and N 2s. The peak at -2.4 ± 0.2 eV arises from the O 2p orbitals. These observations 
are in agreement with that of Zhuravlev and Korabel'nikov [12] that the upper valence bands 
are formed by complex anion states. In the conduction bands, our calculated peaks of 3.25 eV, 
8.25 eV, 9.5 eV, 10.50 eV and 11.25 eV compare favorably with the room temperature vacuum 
ultraviolet (VUV) data from measured values of Yamashita and Kato [37] and Ashida et al. [38], 
respectively. 

A distinct feature of our calculated band structure in Figure 1 is the flat nature of the top of 
the valence band as well as that of the bottom of the conduction band. This singular feature 
means that the band structure of ferro-NaNC>2 can well be described by localized molecular 
levels of an ionic crystal. This observation is similar to that of Ravindran et al. [2], Wang et 
al. [11] and of Jiang et al. [14]. Above the Fermi energy, as can be seen from the DOS, many 
conduction bands exist and some are highly dispersed. From the pDOS for the upper valence 
and the conduction bands, we can infer that the Na states are high up in the conduction bands 
and do not partake in the formation of the molecular solid. 

The top- most valence bands and bottom- most conduction bands are flat. Consequently, the 
exact maximum and minimum values of the valence and conduction bands, respectively, cannot 
be easily identified by merely looking at the plots. From our outputs, before setting the top of 
the valence bands equal to zero, we found that the maximum of the valence band (VBM) is at 
the W point while the minimum of the conduction band (CBM) is at R, respectively, giving a 
fundamental, indirect energy band gap of 2.83 eV from Calculation IV. Our calculated direct 
gaps are 2.90, 2.98, 3.02, 3.22, and 3.51 eV at R, W, X, T, and T, respectively. The scanning 
electron microscopy measurements of Balabinskaya et al. [9] reported a direct band gap of 3.0 
eV. 

The extent to which our electronic structure calculations take temperature into account is 
practically limited to the use of lattice constants obtained at a given temperature. With the 
lattice constants obtained at 120 K, our calculated direct gap of 3.22 eV compares favorably 
with the infrared absorption result of 3.22 eV obtained by Sidman [6] at a low temperature (77 
K). Predictably, it is larger than the room temperature (293.15 K) measurement of 3.14 eV of 
Verkhovskaya and Sonin [7]. 

The positions of our VBM (at W) and CBM (at R) are the same as those from the work 
of Ravindran et al. [11] that found a calculated, indirect band gap of 2.2 eV. Jiang et al. [14] 
calculated an indirect gap of 2.95 eV. However, these authors obtained their VBM at S. In our 
calculations, only the Is states of Na, N, and O, respectively, were in the core. Jiang et al. [14] 
placed Is, 2s, and 2p electrons of sodium in the core. 

Our calculated electron effective masses of 1.18, 0.63, and 0.73 m in the T-X, T-R, and T- 
W directions, respectively, show the highly anisotropic nature of this material. We know of 
no experimental or calculated electron effective masses for NaN02, except the semi-empirical, 
temperature dependent result of El-Dib and Hassan [8] who found an electron effective mass of 
0.61 m Q in the T-R direction. Our result of 0.63, from T-R, agrees with their finding. 
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V. CONCLUSION 



The above results and related discussions are expected to add to our understanding of ferro- 
NaNC>2. The noted agreement between our calculated results and experiment will hopefully 
motivate further experimental and theoretical studies of this material. Our calculated eigenvalues 
at high symmetry points should aid in some comparison with future studies. 

This work and similar ones from this group [17-24] indicate that LDA can correctly describe 
and predict electronic and related properties of semiconductors, provided one methodically search 
for the optimal basis set that minimizes the occupied energies. The predictive capability of such 
calculations is expected to inform and to guide the design and fabrication of semiconductor 
based devices. For example, for binary to quaternary systems, theoretical studies of variations 
of the concentration of an element (or more) have the potential to impact positively industrial 
activities. 
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